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I. INTRODUCTION 



Since its the discovery of interlayer exchange couphng [|I| , this phenomenon has stim- 
ulated a number theoretical investigations. Various approaches have used. These are: 

• the Ruderman-Kittel-Kasuya-Yosida (RKKY) model 

• the free-electron model 

• the hole confinement model |T^|TT[]: 



the Anderson (or srf-mixing) model |12-14 



ab initio calculations [15-30 



The mechanism which now widely accepted for the lEC is based upon quantum in- 
tereferences in the spacer layer due to spin-dependent confinement [p!0| , pTl -|33[| . 

These lecture notes are organized as follows. In, Section I, the quantum inferences 
due to confinement are discussed, and it is shown how this yields an oscillatory interlayer 
exchange coupling. In Section II, the theoretical basis of the one-electron approach of 
Section I is presented. Finally, in Section III, the effect of substitutional disorder is 
addressed. 



II. PHYSICAL MECHANISM OF INTERLAYER EXCHANGE COUPLING IN 
TERMS OF QUANTUM INTERFERENCES: A HEURISTIC APROACH 

The purpose of this section is to present as simply as possible the mechanism of 
interlayer exchange coupling in terms of quantum interferences due to electron confinement 
in the spacer layer. The emphasis here will be on physical concepts rather than on 
mathematical rigor. This discussion is based on the one given in Ref. P2|. 



A. Elementary discussion of quantum confinement 



For the sake of clarity, we shall first consider an extremely simplified model, namely the 
one-dimensional quantum well, which nevertheless contains the essential physics involved 
in the problem. Then, we shall progressively refine the model in order to make it more 
realistic. 

The model consists in a one- dimensional quantum well representing the spacer layer (of 
potential = and width D), sandwiched between two "barriers" A and B of respective 
widths La and Lb, and respective potentials Va and Vb- Note that we use the term 
"barrier" in a general sense, i.e., Va and Vb are not necessarily positive. Furthermore, 
the barrier widths. La and Lb, can be infinite. 



Change of the density of states due to quantum interferences 



Let us consider an electron of wavevector A;+ (with A;+ > 0) propagating to the right 
in the spacer layer; as this electrons arrives on barrier 5, it is partially reflected to the 
left, with a (complex) anplitude tb = \rB\e'^^ ■ The reflected wave of wavevector k~ is 
in turn reflected on barrier A with an amplitude ta = {taI^^'^^, an so on.Q The module 
of the reflection coefficient expresses the magnitude of the reflected wave, whereas 
the argument 4>a{b) represents the phase shift due to the reflection (note that the latter 
is not absolutely determined and depends on the choice of the coodinate origin). 

The interferences between the waves due to the multiple reflections on the barriers 
induce a modiflcation of the density of states in the spacer layer, for the electronic state 
under consideration. The phase shift resulting from a complete round trip in the spacer 
is 

A^ = qD + (PA + (pB , (1) 

with 

q = k+ -k- . (2) 

If the interferences are constructive, i.e., if 

A0 = 2mT (3) 

with n an integer, one has an increase of the density of states; conversely, if the interfer- 
ences are destructive, i.e., if 

= (2n + 1)71 (4) 

one has a reduction of the density of states. Thus, in a flrst approximation, we expect 
the modiflcation of the density of states in the spacer, An{6), to vary with D like 

cos {qD + (j)A + 4>b) ■ (5) 

Furthermore, we expect that this effect will proportional to the amplitude of the reflections 
on barriers A and B, i.e., to Ir^rsl; flnally, An{e) must be proportional to the width D 
of the spacer and to the density of states per unit energy and unit width. 

Tide 



^ Of course, for the one-dimensional model, one has k = — A;^; however, this property will 
generally not hold for three-dimensional systems to be studied below. 



which includes a factor of 2 for spin degenaracy. We can also include the effect of higher 
order interferences due to n round trips in the spacer; the phase shift A0 is then multiplied 
by n and jr^r^l is replaced by jr^r^l". Gathering all the terms, we get, 

2D dg °° 

/^n{e) ^ -tY^ kA'^sT cos n{qD + 0a + (ps) 

(7) 

TT \ de l — rATBC^'^ ) 

As will appear clealy below, it is more convenient to consider the integrated density of 
states 

N{e) = r n{e') de'. (8) 

J —oo 

The modification AN{e) of the intergated density of states due to electron confinement 



IS 



AN{e) = - Im V ^^^^^ e™''^ 

^ n=l ^ 

2 



Im In f 1 - r^r^ e'^^) (9) 

A simple graphical interpretation of the above expression can be obtained by noting that 
Im ln(2;) = Arg (z), for z complex; thus, AN{e) is given by the argument, in complex 
plane, of a point located at an angle A0 = qD + 0^ + 0b on a circle of radius IrA^sl 
centred in 1. This graphical construction is shown in ^ 

The variation of AN{e) as a function of D is shown in ^, for various values of the 
confinement strength IrAr^l- For weak confinement (a), AN{e) varies with D in sinusoidal 
manner. As one increases the confinement strength (b), the oscillations are distorded, 
due to higher order interferences. Finally, for full confinement (c), AN{e) exhibits some 
jumps that correspond to the bound states. We note however, that the period A of 
the oscillations of AN{e) does not depend on the confinement strength, but only on the 
wavevector q = — , namely, A = 27r/g. 

So far, we have implicitely restricted ourselves to positive energy states. Negative 
energy states (i.e., of imaginary wavevector) are forbidden in absence of the barriers A 
and B, because their amplitude diverges either on the right hand side or on the left 
hand side, so that they cannot be normalized. This matter of fact no longer holds in the 
presence of the barriers if Va (or Vb or both) is negative: the negative energy states, i.e., 
varying exponentially in the spacer, can be connected to allowed states of A or B. In 
order to treat these states consistently, we simply have to extend the concept of reflection 
coefficient to to states of imaginary wavevector, which is straightforward. One can check 
that, with this generalization, ^ acounts properly for the contribution of the evanescent 
states. Physically, this can be interpretated as a coupling of A and B by tunnel effect. 
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FIG. 1. Graphical interpretation of equation (^). 
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FIG. 2. Variation of AA^(e) as a function of D, for various values of the confinement strength: 
(a) Iry^r^l = 0.1, (b) Iryir^l = 0.8, (c) jr^r^l = 1 (full confinement). Note the different scales 
along the ordinate axis. 



Energy associated with the quantum interferences in the spacer 



Let us now study the modification of tfie energy of the system due to the quantum 
interferences. In order to conserve the total number of electron, it is convenient to work 
within the grand-canonical ensemble, and to consider the thermodynamic grand-potential, 
which is given by 



1 -I- exp 
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(10) 



At T = 0, this reduces to 



$ = ( ^ {e-ef) n{e) 

J —oo 

= - /'^ N{e) de. 

J —oo 



de 



(11) 



The energy A£' associated with the interferences is the contribution to $ corresponding 

to A7V(£), 



A£; = - Im 

TT 



In (l 



taTb e 
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Three-dimensional layered system 

The generalization of the above discussion to the more realistic case of a three- 
dimensional layered system is immediate. Since the system is invariant by translation 
parallely to the plane, so that the in-plane wavevector ky is a good quantum number. 
Thus, for a given kj|, one has an effective one-dimensional problem analogous to the one 
discussed above. The resulting effect of quantum intereferences is obtained by summing 
over k|| . The modification of integrated density of states per unit area is 

AA^(£) ^-^ImJ d^ky In (l - rArse^^^^) , (13) 

and the coupling energy per unit area is 

AE = ^lmJ d^ky y^°° f{e) In (l - taTb e'^^^) de . (14) 



Quantum size effect in an overlayer 

The case of a thin overlayer deposited on a substrate is of considerable interest. In 
this case, one of the barriers (say. A) consists of the vacuum, and barrier B is constituted 
by the substrate itself. The potential of the vacuum barrier is Kac — ep + W, where W is 



the the work function; thus it is perfectly reflecting for occupied states, i.e., |rvac + 1|- On 
the other hand, the reflection on the substrate (or coefficent rsub) may be total or partial, 
depending on the state under consideration. 

The spectral density of the occupied states in the overlayer can be investigated ex- 
perimentally by photoemission spectroscopy; in addition, by using inverse photoemission, 
one can study the unoccupied states. If furthermore these techniques are used in the 
"angle-resolved" mode, they give information on the spectral density locally in the ky 
plane. 

For a given thickness of the overlayer, the photoemission spectra (either direct or 
inverse) exhibit some maxima and minima corresponding, respectively, to the energies for 
which the interferences are constructive and destructive. When the confinement is total, 
narrow peaks can be observed, which correspond to the quantized confined states in the 
overlayer, as was pointed out by Loly and Pendry [p4 |. 

Quantum size effects due to electron confinement in the photoemission spectra of 
overlayers have been observed in various non-magnetic systems PB|-^ In particular, the 
systems Au(lll)/Ag/vacuum and Cu(lll)/Ag/vacuum offer excellent examples of this 
phenomenon [^ , ^ . 



Paramagnetic overlayer on a ferromagnetic substrate: Spin-polarized quantum size effect 

So far our discussion concerned exclusively non-magnetic systems. Qualitatively new 
behavior can be expected some of the layers are ferromagnetic, case of particular interest 
is the one of a paramagnetic overlayer on a ferromagnetic substrate. 

In the interior of the overlayer, the potential is independent of the spin; therefore the 
propagation of electrons is described by a wave vector fc^ which is spin- independent. The 
reflection coefficient on the vacuum barrier, rvac, is also spin-independent. However, the 
ferromagnetic substrate constitutes a spin-dependent potential barrier; thus, the substrate 
reflection coefficients for electrons with a spin parallel to the majority and minority spin 

T I 

directions of the substrate, respectively r^^^ and r^^j^. It is convenient to define the spin 
average 



rsub = ""^"'"'.^ ""^'"^ (15) 



and the spin asymmetry 

T J. 

Ar,,b = (16) 

In this case, the electron confinement in the overlayer gives rise to a spin-dependent 
modulation of the spectral density versus overlayer thickness; the period of the modulation 
is the same for both spins, whereas the amplitude and phase are expected to be spin- 
dependent. 



The quantum size effects in paramagnetic overlayers on a ferromagnetic substrate have 
been investigated by several groups @-|55|. The systems studied most are Cu overlayers 
on a Co(OOl) substrate and Ag overlayers on a Fe(OOl) substrate. Ortega and Himpsel 
P5| , |i6[] observed a quantum size effect in the normal-emission photoelectron spectra of 
copper overlayer on fee cobalt (001) substrate. They observed peaks due to quantum size 
effects both in the photoemission and in the inverse photoemission spectra an oscillation 
of the photoemission intensity. These quantum size effects manifest themselves also by 
an oscillatory behavior of the photoemission intensity at the Fermi level; as the observed 
oscillation period (5.9 atomic layers) is close to the long period of interlayer exchange 
coupling oscillations in Co/Cu(001)/Co, they pointed out that the two phenomena are 
related to each other; they also claimed that the observed oscillations in photoemission are 
spin dependent and due mostly to minority electrons. A direct confirmation of this con- 



jecture has been given independently by Garrison et al. |^8| and by Carbone et al. |^ by 
means of spin-polarized photoemission. They found that both the intensity and the spin- 
polarization exhibit an oscillatory behavior with the same period (5-6 atomic layers), 
but they have opposit phases, which indicates that the quantum size effect does indeed 
take place predominantly in the minority-spin band as proposed by Ortega and Himpsel 
p5| , |i6[] . Recently, Klasges et al. |Q have observed spin-polarized quantum size effects 
in a copper overlayer on cobalt (001) for a non-zero in-plane wave vector corresponding 
to the short period oscillation of interlayer exchange coupling in Co/Cu(001)/Co; they 
observed short period oscillations of the photoemission intensity in good agrement with 
the short period oscillations of interlayer coupling. This observation provides a further 
confirmation of the relation between quantum size effects in photoemission and oscillation 
of interlayer exchange coupling. 

Photoemission studies of quantum size effects have also been performed in other kinds 
of systems such as ferromagnetic overlayer on a non-magnetic substrate, or systems com- 
prising more layers p6|-p0[ . 

Photoemission spectroscopy undoubtly constitutes a method of choice for investigating 
quantum size effects in metallic overlayers: this is due to its unique features, which allow 
selectivity in energy, in-plane wave vector, and spin. 

Besides photemission, spin-polarized quantum size effects in paramagnetic overlayers 
on a ferromagnetic substrate are also responsible for oscillatory behavior (versus overlayer 
thickness) of spin-polarized secondary electron emission pT| , |62| , linear p3|-p8| and non- 



linear |6^J7^ magneto-optical Kerr, and magnetic anisotropy [n.,72]. However, these 
effects usually involve a summation over all electronic states, so that the quantitative 
analysis of the quantum size effects may be fairly complicated. 



B. Interlayer exchange coupling due to quantum interferences 



Let us now consider the case of a paramagnetic layer sandwiched between two ferro- 
magnetic barriers A and B. Now, the reflection coefficients on both sides of the param- 
agnetic spacer layer are spin dependent. A priori the angle 6 between the magnetizations 
of the two ferromagnetic can take any value; however, for the sake of simplicity, we shall 
restrict ourselves here to the ferromagnetic (F) configuration (ie., 6 = 0) and the antifer- 
romagnetic (AF) one (i.e., 6 = tt). 

For the ferromagnetic configuration, the energy per unit area due to quantum inter- 
ference is easily obtained from 14, i.e., 
1 



Im / d^kn / f{e) In (l - rVLe''^^ ) + In ( 1 - rV^e'^^^ 



de . (17) 



In this equation, the first and the second term correspond respectively to majority- and 
minority-spin electrons. The antiferromagnetic conguration is obtained by reversing the 
magnetization of B, i.e., by interchanging and r^; thus the corresponding energy per 
unit area is 

1 



AE 



AF 



47]-3 
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fie) [in (l 
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Thus, the interlayer exchange coupling energy is 



Ef — Eaf 



1 



Im / d^kii 



f{e) In 
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(1 
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which can be simplified as 

1 f 

Ef - Eaf ^ - ^ Ira d\i / f{e) Ar^Ar^ e'"^^ de 
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de (19) 



(20) 



in the limit of weak confinement. The above expression for the lEC has a rather transpar- 
ent physical interpretation. First, as the integrations on ky over the first two-dimensional 
Brillouin zone and on the energy up to the Fermi level show, the lEC is a sum of contri- 
butions from all occupied electronic states. The contribution of a given electronic state, 
of energy e and in-plane wavevector ky, consists of the product of three factors: the two 
factors ArA and Ar^ express the spin- asymmetry of the confinement due to the magnetic 
layers A and B, respectively, while the exponential factor e*''^^ describes the propagation 
through the spacer and is responsible for the interference (or quantum size) effect. Thus, 
this approach establishes an explicit and direct link between oscillatory lEC and quantum 
size effects such as observed in photoemission. 



Asymptotic behavior for large spacer thicknesses 

In the limit of large spacer thickness D, the exponential factor oscillates rapidly with 
e and ky , which leads to some cancellation of the contributions to the lEC due to the dif- 
ferent electronic states. However, because the integration over energy is abruptly stopped 



at Efi states located at the Fermi level give predominant contributions. Thus the integral 
on e may be calculated by fixing all other factors to their value at ep-, and by developing 
= /cjj^ — A; J around Sp, i.e., 



g±~g±F + 2^-^, (21) 



with 

2 1 1 



[22) 



'^±F ^±F ^±F 

The integration (see Ref. |^ for details) yields 



xF{27ikBT D/hv^p), (23) 



where 



In the above equations, q±F is a vector spanning the complex Fermi surface; the velocity 
v~\^p is a combination of the group velocities at the extremities k~\^p and /cj^. 

Next, the integration on k|[ is performed by noting, that, for large spacer thickness D, 
the only significant contributions arise from the neighboring of critical vectors kj|' where 
q±p is stationary. Around such vectors, q±p may be expanded as 

where the cross terms have been canceled by a proper choice of the and Ky are 

combinations of the curvature radii of the Fermi surface at (k", fc^") and (k", fcj°). 

The integral is calculated by using the stationary phase approximation, and one 
obtains 

xF{27ikBTD/hvl), (26) 
where g", f", Ar^, Arg correspond to the critical vector k|, and 

)^/^(-^)'''; (27) 



in the above equation, one takes the square root with an argument between and vr. 

This analysis shows that in fine, the only remaining terms in the limit of large spacer 
thickness D arise from the neighborhood of states having in-plane wavevectors kjj' such 
that the spanning vector of the Fermi surface q±F = k^p ~ ^If is stationary with respect 



FIG. 3. Sketch showing the wavevector giving the oscillation period of oscillatory inter- 
layer exchange coupling, for the case of a non-spherical Fermi surface. 



to k|| for k|| =k°, and the corresponding contribution oscillates with a wavevector equal 
to q'^p. This selection rule was first derived in the context of the RKKY model 0; it is 
illustrated in Fig. ^. There may be several such stationary spanning vectors and, hence, 
several oscillatory components; they are labelled by the index a. 

The above selection rule allows to predict the oscillation period(s) of the interlayer 
exchange coupling versus spacer thickness by just inspecting the bulk Fermi surface of the 
spacer material. In view of an experimental test of these predictions, noble metal spacer 
layers appear to be the best suited candidates; there are several reasons for this choice: 

• Fermi surfaces of noble metals are known very accurately from de Haas-van Alphen 
and cyclotron resonance experiments 



• since only the sp band intersect the Fermi level, the Fermi surface is rather simple, 
and does not depart very much from a free-electron Fermi sphere; 

• samples of very good quality with noble metals as a spacer layer could be prepared. 

Fig. ^ shows a cross-section of the Fermi surface of Cu, indicating the stationary 
spanning vectors for the (001), (111), and (110) crystalline orientations [^]; the Fermi 
surfaces of Ag and Au are qualitatively similar. For the (111) orientation, a single (long) 
period is predicted; for the (001) orientation, both a long period and a short period 
are predicted; for the (110) orientation, four different periods are predicted (only one 
stationary spanning vector is seen in figure ^, the three others being located in other 
cross-sections of the Fermi surface). These theoretical predictions have been confirmed 
successfully by numerous experimental observations. In particular, the coexistence of a 
long and a short period for the (001) orientation has been confirmed for Cu ||7^J75[] , Ag 




FIG. 4. Cross section of the Fermi surface of Cu along the (110) plane passing through 
the origin. The solid dots indicate the reciprocal alttice vectors. The dashed lines indicate 
the boundary of the first brilouin zone. The solid arrows, respectiveley horizontal, oblique, and 
vertical, indicate the vectors giving the oscillation period(s), respectively for the (001), (111), 
and (110) orientations. 



and Au [[T^TDI; and the experimental periods have been found to be in excellent 
agreement with the theoretical ones. 

In a further attempt to test the theoretical predictions for the periods of oscillatory 
coupling, several groups [p^-p6[] have undertaken to modify in a controlled manner the size 



of the Fermi surface (and hence, the period of the coupling) by alloying the spacer noble 
metal (Cu) with a metal of lower valence (Ni); in both cases, the change in oscillation 
period due to alloying has been found in good agreement with the expected change in the 
Fermi surface. 

Although the asymptotic approximation is often an excellent one to describe the lEC, 
in cases where the reflection coefficients vary strongly with e and ky near the station- 
ary spanning vectors of the spacer Fermi surface, preasymptotic corrections nead to 
be considered. A detailled discussion of the preasymptotic corrections is presented in 
Ref. 



Effect of magnetic layer thickness 

As already mentioned, the influence of the lEC on the ferromagnetic layer thickness 
is contained in the reflection coefficients Ar^ and Ar^. If the ferromagnetic layers are 
of finite thickness, reflections usually may take place at the two interfaces bounding the 



TABLE L Comparison between the theoretical predictions of Ref. [^and experimental ob- 
servations for the oscillation periods of interlayer exchange coupling versus overlayer thickness. 
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FIG. 5. Sketch of the waves 
layer of finite thickness L. 



contributing to the net reflection coefficient on a ferromagnetic 
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FIG. 6. Contour plot of the interlayer exchange couphng constant J = [Ep — Eaf)/"^ vs 
spacer thickness D and magnetic layer thickness L, calculated within the free-electron model 
(see Ref. [32| for details). The spacing between successive contour lines is 40 x 10^'^ erg cm~'^; 
the shaded area corresponds to antiferromagnetic coupling. 



ferromagnetic layers, as sketched in Fig. giving rise to interferences and hence, to 
oscillations of the lEC versus ferromagnetic layers thickness. A more detailled discussion 



of this effect is given in Refs. |3^j8^. This behavior was first predicted from calcula- 
tions based upon a free-electron model 0. On the experimental point of view, it was 
confirmed by Bloemen et al. in Co/Cu/Co(001) and by Okuno and Inomata in 
Fe/Cr/Fe(001). The amplitude of the oscillations of the lEC versus ferromagnetic layers 
thickness is generally much smaller than the oscillations versus spacer thickness, and do 
not give rise to changes of sign of the lEC. This is illustrated in Fig. 0. 



Effect of overlayer thickness 

A more (a priori) surprising behavior is the dependence of the lEC on the thickness 
of the protective overlayer. From a naive point of view, one might think that layers 
external to the basic ferromagnet /spacer /ferromagnet sandwich should not influence the 
interaction between the two ferromagnetic layers. This view is incorrect, in particular 
when the system is covered by an ultrathin protective overlayer. In this case, the electrons 
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FIG. 7. Sketch of the model used to discuss the influence of the overlayer and vacuum barrier 
on the reflection coefficient. 

are able to reach the vacuum barrier, which is a perfectly reflecting one, so that strong 
confinement and interference effects take place in the overlayer, which lead to a weak but 
sizeable oscillatory variation of the lEC as a function of the overlayer thickness. 

Here, we illustrate the effect of interferences in the overlayer by means of the simple 
model depicted in Fig. ^ An electron of energy e and in-plane wave vector bf ky = 
is incoming from the right. The magnetic layer is represented by a barrier of height V 
and width L; the overlayer has a thickness T, and the vacuum is modelized by a semi- 
inifinite potential barrier of height W . The spin dependence of the magnetic barrier is 
not considered here. 

Since the energy of the incoming electron is smaller than the vacuum barrier height, 
particle flux conservation imposes to have |r| = 1, and the reflection coefficient may be 
written as r = e^"^. It is thus sufficient to discuss the behavior of the reflection phase 
shift (f). The variation of as a function of the overlayer thickness T is shown in Fig. 
The behavior is very contrasted, depending on the value of V . Two limit cases can be 
considered. 

In the first case, V is small as compared to e, the reflection phase shift is essentially 
given by the sum of the phase shift due to the reflection on vacuum and and of the one 
due to the round trip through the magnetic layer and the overlayer, so that one has a 
linear variation of versus T, as appears in Fig. ^ for V = Q and V = 0.5 eV. Thus, the 
reflection coefficient r = e**^ varies periodically versus overlayer thickness T, with a period 
equal to 2T[/q±. The change of phase shift as V varies is due to the change in the phase 
shift associated with the travel through the magnetic layer. 

The opposit limit case is obtained when e <V\ m this case, r is essentially constant. 




FIG. 8. Variation of the reflection phase shift (j) versus overlayer thickness T for various values 
of the height V of the magnetic barrier, as indicated by the arrows. Parameters: £ = 1.5 eV, 
W = 5.0 eV, L = 6.0 A. 
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FIG. 9. Sketch of the model used to discuss the influence of the over layer and vacuum 
barrier on the inter layer exchange coupling; the solid and dashed line represent, respectively, 
the majority spin and minority spin potential profiles, for the configuration of ferromagnetic 
alignement. 



exepct for resonances where (j) makes a jump of 27r, corresponding to the crossing quasi- 
bound states in the overlayer. This clearly demonstrated in Fig. for V = 2.0 eV and 
V = 2.5 eV. As the transmission through the magnetic barrier decreases, the resonances 
become narrower and their effect is completely negligible. 

The intermediate situation evolves continuously between the two limit the 
curves corresponding to V = 1.0 eV and = 1.5 eV in Fig. |^ show. 

With help of the above argument, we can understand easily how the thickness of 
the overlayer influences the lEC As will be discussed below, the variation of the inter- 
layer coupling versus overlayer thickness T has a completely different behavior, depending 
on wether the magnetic barrier strongly confines the minority spin electrons at Fermi 
level(i.e., ep < A^x) or not; for brievety, in the following these two opposit situations will 
be refered to as strong confinement and weak confinement, respectively. 

Let us consider first the case of weak confinement. In this case, one finds that the lEC 
dependence upon overlayer thickness T is via an oscillatory function ofT + D, as appears 
clearly from Fig. |10| (top). 

On the other hand, for the strong confinement Fig. (bottom) shows, the 

interferences in the overlayer have only a small influence and yields a small modulation 
of the coupling strength as T varies. 

This effect, which follows directly from the quantum interference (or quantum size 
effect) mechanism, has been proposed and experimentally confirmed independently by de 
Vries et al. |^T| for the Co/Cu/Co(001) system with a Cu(OOl) overlayer, by Okuno and 
Inomata for the Fe/Au/Fe(001) system with a Au(OOl) overlayer, and by Bounouh 
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FIG. 10. Contour plot of the interlayer exchange couphng, Ep—Ep^Yi versus spacer thickness 
D and overlayer thickness T. Parameters: top pannel: ep = 1-5 eV, W = 5.0 eV, Aex = 0.5 eV, 
L = 6.0 A; bottom pannel: ep = 1-5 eV, W = 5.0 eV, Aex = 2.0 eV, L = 6.0 A. The spacing 
between succesive contour lines is repectively 2.0 x 10~^ erg.cm"^ and 0.2 erg.cm"^ in the top 
and bottom pannels; the shaded area corresponds to antiferromagnetic coupling. 



TABLE IL Comparison between the theoretical predictions of Ref. and experimental 
observations for the oscillation periods of interlayer exchange coupling versus overlayer thickness. 



overlayer 


theoretical periods 


system 


experimental periods Ref. 


Cu(OOl) 


Ai = 2.6 AL 


Cu/Co/Cu/Co/Cu(001) 


A 


5 AL 




A2 = 5.9 AL 








Au(OOl) 


Ai = 2.5 AL 


Au/Fe/Au/Fe/Au(001) 


Ai « 


2.6 AL II 




A2 = 8.6 AL 




A2 « 


8.0 AL 


Au(lll) 


A= 4.8 AL 


Au/Co/Au/Co/Au(lll) 


A ^ 


5 AL II 



et ai [0 for the Co/Au/Co(0001) with a Au(lll) overlayer. In both cases, the observed 
period(s) for the oscillations versus overlayer thickness were found to be in good agreement 
with the theoretically predicted ones. This effect has also been confirmed by means of first- 
principles calculations for the Co/Cu/Co(001) with a Cu overlayer The comparison 
between the periods of oscillations versus overlayer thickness predicted theoretically and 
those observed experimentally is given in Table 



III. EXCHANGE INTERACTIONS WITHIN DENSITY FUNCTIONAL 

THEORY 

The mechanism of lEC presented above is based upon an independent electron pic- 
ture. That such a picture is actually valid is not immediatly obvious, in view of the fact 
that exchange interactions are indeed due to the Coulomb interaction between electrons. 
In this Section, we present the justification of the independent electron picture for the 
lEC It relies on the use of "force theorems" for the magnetic interactions, which follow 
from variational properties of various energy functionals. Since these considerations are 
not particular to the problem of interlayer coupling, but apply to any kind of exchange 
interactions in metallic magnetic systems, a more general point of view will be adopted 
here. 



A. Constrained density functional theory 



Almost all modern methods of electronic structure calculation rely the density func- 
tional theory (DFT) of Hohenberg and Kohn |Q , within the local density approximation 



(LDA) of Kohn and Sham ||9^. In the DFT (generalized to take the spin polarization into 
account |9^) the system under consideration is described in terms of the local density 
spinor n(r) whose matrix elements are given byQ 



^Throughout this paper the tilde ~ will be used to note spinor quantities. 




In the above equation ipl^{r) (^/'^(r)) is the second-quantization creation (destruction) 
operator for an electron of spin a {j3) at r, and is the many-body wave function. 
Alternatively, the system can be described in terms of the more familiar local charge and 
spin densities, p{r) and m(r), which are related to the density spinor by 

p = tr (n) (29) 
m = tr (n or) (30) 
n = (P^^+^^-^^ (31) 

where 
and 

(33) 

is the vector whose components are the Pauli matrices; x, y and z are the unit vectors of 
the cartesian axes, and "tr" represents the trace over spin indices. 

The DFT establishes that the total energy of a given system is a unique functional 
S[n\ of the density spinor n{r) which is stationary and has its absolute minimum for 
the density spinor fio(r) corresponding to the ground state of the system, and that the 
minimum energy is equal to the ground state energy, £q, i.e., 

So = S [no] = min nS[n]- (34) 

The approach universally used to compute the ground state density spinor no and the 
ground state energy Sq is the one proposed by Kohn and Sham |9^, who rewrote the 
energy functional as 

S [n] = To [n] + Vext [n] + Kh [n] + ^xc [n] • (35) 

In the above equation, the first term is the kinetic energy of a fictitious system of inde- 
pendent electrons having the same density spinor n(r) as the real system. The second 
and third terms are respectively the potential energy of the electrons in the external 
potential and the Hartree approximation to the energy of Coulomb repulsion between 
electrons. The last term represents the exchange and correlation correction to the kinetic 
and Coulomb energies. 

The Hartree and external potential terms can be calculated trivially. The exchange- 
correlation term cannot be calculated exactly; however, convenient and efficient approxi- 
mations exist for computing it within the LDA. The most problematic term is the kinetic 



energy 7o[n], for which no satisfactory approximation is known. Kohn and Sham solved 
the problem by showing that the system can be mapped to a fictitious system of indepen- 
dent electrons having the same density spinor and moving in an effective local potential 
spinor given by 

WeS [n] (r) = Wext (r) + Wh [n] (r) + w^c [n] (r) (36) 
where tyextlr) is the external potential spinor, and 

" w ^''^ 

are respectively the Hartree and exchange-correlation potential spinors. One has to solve 
Schrodinger-like equations for independent electrons in the effective potential spinor Wefr[^] 
(Kohn-Sham equations). The ground state energy is then given by 

£0= J2 + "^'ext [^0] + l^H [no] + £xc [no] - ^d^r tr (no Wefr [^0] ) (39) 

£i<£F 

where Ei are the single particle energies corresponding to the solutions of the Kohn-Sham 
equations. Since the effective potential spinor depends on the density spinor, the Kohn- 
Sham equations must be solved self-consistently, which is usually achieved iteratively, 
starting from a trial density spinor. 

As it stands, this theory is not very convenient for computing exchange interactions: 
indeed, as it gives some information only on the ground state of a system, one would learn 
only which configuration of magnetic moments (e.g., ferromagnetic, antiferromagnetic, or 
canted) corresponds to the ground state; thus would obtain the sign of the interaction, but 
one would not learn anything about its strength. For this, one would have to include in 
the Hamiltonian the Zeeman interaction due to an external magnetic field, and compute 
the configuration of magnetic moments, and the total energy of the system as a function 
of this external field. This precisely how one proceeds in an experiment! Although this is 
a conceptually straightforward approach, it is not very convenient to implement, and has 
never been actually used. 

Actually, what we would like to know is the ground state energy of the system, sub- 
jected to the restriction that the local spin-polarization m(r) is constrained to be along 
some prescribed direction of unit vector u(r); the constraint may be extended to the 
whole space or restricted to a given subspace V. This approach relies on a particular case 
of the constrained density functional theory (CDFT) of Dederichs et al [^]. It uses the 
standard method of Lagrange transformation, in which one defines the new functional 

^u[n, hxu] =m - I d'r m(r) ■ h^fl(r). (40) 

The Lagrange parameter is a transverse magnetic field, perpendicular to the local 
magnetization direction u, to be determined self-constently. The density spinor of the 



constrained ground state n^(r, its energy fexch[u], and the corresponding transverse field 
are obtained by minimizing J^aln, h^a] with respect to n and h^a, i-e., 

S[nl] = min^i^h^fl ^u[^, h^a] = ^u[^a> ^la] (41) 

= ^exch[u] (42) 

The physical meaning of hxa(r) is quite clear: it is the transverse external field one needs 
to adjust in order to maintain everywhere the magnetization parallel to the prescribed 
direction u(r). The local density of torque due to exchange interactions upon imposing 
the constraint u(r) is given by 

where the vector r2(r) represents a local rotation of the spin-polarization axis. It is related 
to the transverse field by 

ra(r) = -ma(r) xhxa- (44) 

The constrained ground state density spinor is calculated as in the unconstrained 
DFT, by solving self-consistently the Kohn-Sham equations with the effective potential 
spinor 

Wes[n, h_La](r) = Wext(r) + WH[n]{r) + Wxc[n](r) - h_La(r) • cr. (45) 
The constrained ground state energy is then calculated as 

^^exch[u] = E - /d'r tr {hi W,sK, Kul) + ^eAKl + ^H[nl] + ^AK] 

ei<ep 

- I d^r mt • hla- (46) 

The functional £exch[u], or, equivalently, the torque density ra(r), contains all the infor- 
mation we may wish to know about exchange interactions in the system. 



B. Harris-Foulkes functional and "force theorems" 

The self-consistent method described above constitute a conceptually straightforward 
approach for computing exchange interactions. However, this approach also has severe 
drawbacks. 

• First of all, self-consistent calculations are computationally very demanding, because 
a large number of iterations are usually required to achieve convergence towards the 
self-consistent solution with sufficient accuracy. 



• Second, within the self-consistent approach the exchange energy, which is typically 
of the order of 10~^ to 10~^ eV/atom, is obtained as the difference between the total 
energies for two different configurations, which are of the order of 10'* eV/atom; thus, 
the total energies must be obtained with a tremenduous relative accuracy in order 
to get reliable results. In this case, the results may be plagued by roundoff errors. 

• But the most serious difficulty we face when using the self-consistent approach is 
that it provides us very little physical insight about the mechanism of exchange 
interaction. 

For all the reasons mentioned above, it is desirable to develop a method which is 
computationally convenient and accurate, and at the same provides a clear physical picture 
of the mechanism of exchange interaction. The fact the energy of exchange interactions is 
small as compared to the total energy of the system suggests us that a perturbation-like 
theory might be appropriate. Alternatively, one can exploit the variational properties of 
the density functional to compute approximately the constrained ground state energy in 
a single non-self-consistent shoot. This the approach we shall adopt here. We shall derive 
a force theorem [p.0CI|J101|| that allows to express £^exch[u] (within an unimportant constant) 



in terms of the sum of single-particle energies, calculated non-self-consistently. 
Because of its variational property of the energy fonctional J^u[^;h_|_a] satisfies 

^u[n, h^u] = ^exch[u] + 02{n- , h^a - h^-) (47) 

where O2 {n — h\, h^Q — h^^) is of second order with respect to n — n\ and h_|_a — h^- . 
So, if we have good guesses n and hj_u for the exact n\ and h^^, and if we are able to 
compute J^u[?^, h^u], we get an estimate of £^exch[u] which is accurate to second order in 
the error of our initial guesses. As we shall see just below, the difficulty of this approach 
lies in the calculation of J^yi[n, hxa] for the chosen n and h^fj. 

This is seen as follows. From a trial input density and input transverse field h^^^ 
one gets the effective potential spinor Wes[nin, h™u]i by solving the Kohn-Sham equations 
for this effective potential spinor one gets a set of single particle energies ei[ni-a_, h™j^], from 
which one can in turn compute an output density fiout- From this, we can compute 



£i<£F 

+Vext[^out] +UH[nont] + ^^xc[^out] " / d^r mout ■ h^ji (48) 

We now see where the problem lies: we get an estimate of the energy fonctional corre- 
sponding to the output density spinor, not to the input one; furthermore, the expression 
mixes in a complicated manner and fiout- Since we don't know the input density spinor 
corresponding to an arbitrary output density spinor, this cannot be used to compute the 
energy fonctional for a chosen density spinor. 



This problem was circumvented by Harris ||102|| who defined an auxihary functional as 
follows 

Gil [n-m, Win, i^±u] = J2 ~ / (^in ^in) + Vcxt [^in] + l^H [^in] + £kc [^in] 

£i<ep 

- I ■ hi\ (49) 



where riin is the input density, £:j[^in] ^'^^ the single particle energies calculated with the 
input effective potenti wi^, and h^j^ is the input transverse field. Here we have used 
a generalization of the Harris functional proposed by Foulkes and Haydock |103| . It is 
straighforward to show that 

e?fl[n, w, h^fl] = J^fl[nt , hi-] + 02 (n - n^, w - wl, h^a - h^-) (50) 

where 

= Wcxt + wnlnl] + w^c[K] - ■ (51) 

and O2 {n — n(^, w — vji^, h_Lu — h^^^) is a (non generally positive) error of second order 
with respect n — n^, w — and h^u — h^u- The properties of the Harris functionals have 
been discussed by a number of authors [|104|-|106| who found that it often yields better 



approximate estimations of the ground state energy than the Hohenberg-Kohn functional. 

We now have all the material needed to establish the force theorem for magnetic 
exchange interactions. To this end, we make the following choice for the input values of 
h, w and h±i^■. 

hTu(r) = (52b) 
w;in(r) = Vinir) ao - /lin(r) u(r) • a. (52c) 

The input charge and spin densities are chosen to be independent of u in magnitude, 
with the axis of the spin-polarization along u; thus we neglect the redistribution of charge 
and magnetic moment due to rotating the magnetic moments, which usually constitutes a 
good approximation. The input effective field is chosen to be parallel to u; this is usually 
a good approximation, because the transverse component of the effective field is much 
smaller than its longitudinal component, which is of the order of the exchange splitting. 

For the magnitude of the input effective potential fin and input effective field /lin we 
can specify a little more by ascribing them the value corresponding to the LDA, i.e., 

Vin{r) = i;cxt(r) + VH[pin]{r) + v^^^[pin, min](r) (53a) 
/lin(r) = /i^DA[pin,min](r). (53b) 

However, this choice is not necessary and the stationarity of the Harris functional allows 
more fiexibility in the choice of the input values for the effective potential and effective 
field. 



Inserting the particular choice ( p2a| ,6,c) for the input quantities in the definition of 
the Harris functional (|49|) we find 

Gu[nm, Win, h'°u] = ^*bin, ^inU] - J d^r (pin^^in " "^in/iin) 

£i<£F 

+Vext[Pm] + W//[piJ + £^c[pm, "^in]- (54) 

We find that the only term which depends on the constraint u is the sum of single particle 
energies. Thus, we finally obtain 

^^exch[Ul] - £^exch[u2] = £i[Vm, KAl] - ^^iKn, ^111^2] 

+C2(Pm-Pfl,min-m^,h*-) (55) 

which consitutes the force theorem for the exchange interaction energy. This is an impor- 
tant result, for several reason: 

• The energy of exchange interactions £^cxch[u] can be computed (within an unim- 
portant constant) in a single shoot, without need for lengthy iterations towards 
self- consistently, which makes the calculations considerably faster and easier. 

• Although the exchange energy is ultimately due to Coulomb interactions between 
electrons, it is expressed here as a sum of single-particle energies for indenpendent 
electrons. This is of considerable practical importance, because the sum of single 
particle energies is much smaller than the total energy of the system; thus the relative 
accuracy needed is much less than for self-consistent total energy calculations, and 
the risk of computational error is much smaller. 

• A further remarkable feature of this result is that the exchange energy difference 
depends only on the input potential and input field, and not at all on the input 
charge and spin densities (as long as the latter don't vary much upon rotating the 
moments). This leaves considerable flexibility in setting up practical computational 
schemes, and allows to use suitably parametrized schemes, without jeopardizing 
seriously the accuracy of the results. 

The torque density is given by 

(5£^cxch[u] 



ra(r) 



5n{vy ^^^^ 



From the Hellman-Feynman theorem, we obtain for a given single-particle energy 
corresponding to the wave function |0j[u]). 



Upon performing the rotation fl{r) = f2(r)nn(r), which changes u(r) into u'(r), the 
effective potential WcfT(i") is transformed according to (for simphcity, the subscript "eff" 
will be ommited below) 

t^u(r) ^ Wu'{r) = R{r) w^ir) R-\r) (58) 

where the rotation matrix is given by 

R{r) = e-'"W-* 

= cos(^^) -ino(r)fi(r)-^sin(^^^. (59) 

Thus, 

6H 



and by using the relation 



|r) /ii,(r) u ■ a (r| (60) 



mout(r)= E {M^]\r)d- (r|0,[u]) (61) 

Si<£F 

we finally obtain 

rfl(r) = - /lin(r) u X mout(r). (62) 
The above result constitutes the force theorem for the exchange torque density. 



C. Examples of first-principles calculations of lEC 

The use of the magnetic "force theorems" discussed in the preceding Section allows 
to reduce by several orders of magnitude the computation times required to calculate the 
interlayer exchange coupling. Such first-principles calculations have been performed by a 
number of authors [p!8|-|30|. 



The Fig. |TT| shows the interlayer exchange coupling for the Co/Cu/Co(001) systems, 
for semi-infinite Co layers (top), and 5 AL thick Co layers (bottom) ||107|| . The data have 
been multiplied by N"^ in order to highlight the dependence of the lEC, as predicted 
from the asymptotic approximation. The data corresponding to the two different Co 
thicknesses seem very different at first sight. In order to evidence more clearly the periodic 
oscillatory behavior, one can perform a discrete Fourier analysis of the N dpendence of 
the lEC I^Tj. The corresponding results are shown in Fig. 

One sees very clearly that the oscillations comprise two superimposed oscillatory com- 
ponents, as expected, with periods in good agreement with the ones predicted in terms of 
the spacer Fermi surface On the other hand the amplitude of the oscillations change 
dramatically, as the Co thickness is varied. 
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FIG. 11. Calculated interlayer exchange coupling N'^Ex as a function of spacer thickness A'^ 
(in AL) for the Co/Cu/Co(001) system; top panel: semi-infinite Co layers; bottom panel: 5 AL 
thick Co layers [|lO§. 




FIG. 12. Discrete Fourier transform of the data shown in Fig. left panel: semi-infinite 
Co layers; right panel: 5 AL thick Co layers |107]. 



These results provide a very clear confirmation of the results obtained from the discus- 
sion in terms of quantum interferences in the asymptotic regime, namely that the periods 
of oscillation versus spacer thickness depend only on the spacer material, but the ampli- 
tudes can vary with the thickness (as in the present case) or nature of the ferromagnetic 
layers. 



IV. SUBSTITUTIONAL DISORDER 



The case where (some part of) the system under consideration consists of a substi- 
tutionally disordered alloys is extremely important in practice. This happens of course 
if some layers are by purpose chosen to consist of alloy; but more generally, unavoidable 
interdiffusion takes place at the interfaces, giving rise to a disordered interface region with 
a progressive variation of concentration. 



A. "Vertex cancellation" theorem 



We present here a general discussion of exchange interactions in the presence of substi- 
tutional disorder. The results given here are used in the present paper to study interlayer 
exchange interactions, but they are also applicable for studying exchange interactions 
within a ferromagnet, exchange stiffnesses, spin- wave energies, etc. 



The principal result is the "vertex cancellation theorem" of Bruno et al. ||108|| . We 
present here an alternative, more general, derivation of this result. 

Let us specify the notations used here. Our purpose is to compute the total energy 
(more precisely the thermodynamic grand-potential) as a function of the prescribed local 
direction u(r) of the magnetization. Explicitely, the Hamiltonian operator corresponding 
to a particular configuration u = {u(r)} of the local moments is written as 

H = K + Va. (63) 

The matrix elements of the kinetic energy operator K are 

(r|K|r')=<5(r-r')^;^ao, (64) 
Im dr^ 

and those of the potential operator V 

(r IVfll r') = 5{j - v') [v{j) - h(j) u(r) • a] , (65) 

where w(r) and h{Y) are respectively the local effective potential and the local effective 
field. The corresponding Green's function operator is 



(66) 



Because of the substitutional disorder, v{r) and h{r) have some spacial randomness. 
However, within the domain of apphcabihty of the "force theorem" , we assume that they 
are unchanged upon changing the moment configuration u. 

The central quantity in the theory of disordered alloys is the configuration averaged 
Green's function 



G^{z) = {Gi,{z))^ = [z-K-T^{z)] 



-1 



(67) 



where {■ ■ indicates an average over all possible alloy configurations, and where ^u.{z) 
is the self-energy. The Green's function is given by 



Gu{z) = Gu(z) + Gu(z)Tu(2;)Gfl(z), 
where the t-matrix Tij^{z) is given by 

T^{z) = (Vfl - U{z)) fl - Ga(^) (Vfl - Uiz))r' 



{61 



(69) 



iFTom the definition of the configuration averaged green's function, we get the self- 
consistency condition 



or equivalently. 



(Tu(^)), = 0, 



1 - Gfl(z) (Va - T^{z)) 



-1 



(70) 



(71) 



In practice, the self-energy satisfying the conditions (^^, |7Tp cannot be calculated exactly 
and one has to resort to approximations. The most popular approach is the CPA, in 
which conditions (|70| , [Til) are replaced by on site conditions, for all atomic sites R. 

If one uses the "force theorem" , the thermodynamic grand-potential includes only the 
single-particule energies (Kohn-Sham eigenvalues) and is given by 



+ 00 



defie)N^ie) 



(72) 



where f{e) is the Fermi-Dirac function, and A^u(^) the integrated density of states (av- 
eraged over all possible alloy configurations) for the local moment configuration u. The 
integrated density of states is given by 



Ar-(e) = ilmTr[(lnG(^))jr+^°" 



'cJ^=-oo+iO+ ■ 



(73) 



Let us consider the quantity 



Aa(^) = Tr(lnGfl(z))^ 



(74) 



By some simple algebra, one can show that 



(75) 



where 



X^{z) = Tr (In [1 - Gfi(z) (Vq - Z^z))])^ 



(76) 



is called the vertex correction. This term is usually non negligible and is difficult to 
calculate. We shall show however that its variation upon varying u takes a simple form. 

Let us precise the variation of the Hamiltonian under varying u. If we chose a reference 
configuration Uq for the local moment configuration, the configuration u is obtained by 
performing locally a rotation of vector 



0(r) = l^(r) nn(r). 

Explicitly, one has 

u = (fin • Uo) Afj + cos Q. (fin x u) x fin + sin Q. (n^ x u). 
The potential operator Va corresponding to u is obtained from as 

and the matrix elements of the rotation operator are given by 

(r |Rn| r') = 6{t - r') exp (-^^^(r) • a-) 

'fi(r 



(77) 



(78) 



(79) 



= 5(r - r') 



cos 



(To - ism — — nr2(r) • cr 



(80) 



We shall also make of the relation 

6 (RnBR^^) 
5Q{r) 



|r)a-(r| ; RnBR^.^ 



(81) 



for an arbitrary operator B. 

The derivative of (from now on, we shall omit the dependence on complex energy 
z) with respect to f2 is given by 



5n{r) 



Tr 



Of 



"(5ri(r) 



(82) 



with 



-Tr 



6Gi 



(Va-Za) l-G<i(Va-Za) 



-Tr G,- 



l-Gfi(Va-Za) 



(83) 



The first term on the right-hand side of the above equation is zero, because of condition 
([70|). Next, we split the self-energy in two parts as 



where T.^-^^ varies with f2 like Vu, i.e.. 



Then, we obtain 



6n(r) 



Tr 



^4) 



^5) 



(86) 



By using the permutation invariance of the trace and the conditions (^, 71), we obtain 
finally 



and hence 



50(r) 



_5A^ 
5n(r) 



Tr 



c u 
"5f^(r) 



^7) 



Tr 



(1) 



"(jrifr) 



Thus, the torque density due to the exchange interactions is given by 



5n{T) 



f{e) Im Tr 



Gfl(£ + iO^ 



5n(r) 



de. 



^9) 



This exact result constitutes the "vertex cancellation theorem" for the torque density. 
Its usefulness arises from the fact that the "vertex corrections" have been eliminated. 
A further important feature is that the exact torque is given in terms of the angular 
derivative of T.^^'' only, which is known explicitly because it follows the angular variation 
of Vu; the angular derivative of the remaining part of the self-energy, 51)^ , which cannot 
be calculated easily, has been eliminated. 

In order to compute the difference of thermodynamic grand-potential between two 
local moment configurations ui and U2, we use a theorem due to Ducastelle ||109|| , which 
states that the thermodynamic grand-potential, considered as a functional ^[G, V] of the 
independent variables G and V, is stationary with respect to G when the latter satisfies 
the self-consistency condition (^, [7T| ). This means that a first-order error on Gq gives 



only a second-order error on $u- Writing Gu as 



(90) 



with 



G^l\z) ^ {z - K - T^^\z: 



(91) 



7-(l) 



we take Gj^ as a trial value for computing $0. This can be expected to be a good 
approximation, provided the condition 



m r 



drj 



dr 



(92) 



(where p(r) and m(r) are respectively the electron and spin densitites) is satisfied. Re- 
placing Gu by G^j in Eq. (^9l), we can now integrate over angles, and we get 



^Ul ^U2 



f{e) Im Tr 



lnGl'\£ + iO+)-lnG£^(£ + iO^ 



(93) 



which constitutes the "vertex cancellation theorem" for exchange energies. In the deriva- 
tion of the "vertex cancellation theorem", we have made use of the exact self-consistency 
condition (|7D|, |7I]) ; one can show also that the same result holds if one uses the approxi- 
mate CPA self-consistency condition. 

In the case of interlayer coupling, the condition (|92|) is satisfied even for large rotation 
angles, because df2/dr differs from zero only in a region where m(r) is negligible. This 



was confirmed by explicit numerical calculations in Ref. [108 



B. Numerical studies of lEC in presence of alloy disorder 



The interlayer exchange coupling through a Cui-^^Ni^ spacer layer has been studied 

These authors have observed a decrease of 



experimentally by several groups |[84-86 



oscillation periods with increasing Ni concentration, which they have attributed to the 
shrinking of the Fermi surface. 

First-principles calculations of lEC for 

the Co/Cui_^.Ni^/Co(001), Co/Cui_^.Zn^/Co(001), and Co/Cui_^Au^/Co(001) systems 
based upon the "vertex cancellation" theorem have been performed by Kudrnovsky et al. 
| 110|| . Their results are shown in Fig. |1^. The systematic shift of the oscillation periods 
with impurity concentration appears clearly. The decrease (resp. increase) of oscillation 
periods with increasing concentration of Ni (resp. Zn) is clearly due to decrease (resp. 
increase) number of conduction electrons. Alloying with an isoelectronic metal (Au) on 
the other, yields quasi- const ant oscillation periods. 




FIG. 13. Left panel: composition dependence of the absolute values of a discrete Fourier 
transform of Ex at T=0 K for two Co(OOl) slabs each five monolayers thick separated by an 
fcc-Cui-ajNi^; alloy spacer: (i) Cuo.75Nio.25 (full line), (ii) Cuo.85Nio.15 (dashed line), Cuo.gNio.i 
(dashed-dotted line), and (iv) an ideal Cu spacer (dotted line); from Ref. [ |110| |. Right panel: 
composition dependence of the coupling periods for two Co(OOl) slabs each five AL thick sepa- 
rated by an fcc-Cui-^M^ alloy spacer: (i) M=Ni (bullets), (ii) M=Au (squares), and (iii) M=Zn 
(diamonds). The lines serves as a guide for the eye and distinguish between short (full lines) and 
long (dotted lines) period oscillations. Open circles for Cuq. 75^0.25 and CU0.5AU0.5 represent 



the approximate virtual-crystal values. The periods are given in ALs; from Ref. 110 |. 
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